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Type la supernovae (SNe la) have been used as excellent standardizable can- 
dles for measuring cosmic expansion, but their progenitors are still elusive. 
Here we report that the spectral diversity of SNe la is tied to their birthplace 
environments. We find that those with high- velocity ejecta are substantially 
more concentrated in the inner and brighter regions of their host galaxies than 
are normal-velocity SNe la. Furthermore, the former tend to inhabit larger 
and more-luminous hosts. These results suggest that high-velocity SNe la 
likely originate from relatively younger and more metal-rich progenitors than 
normal- velocity SNe la, and are restricted to galaxies with substantial chemi- 
cal evolution. 

Type la supernovae (SNe la) are among the most energetic and relatively uniform stellar 
explosions in the Universe, and were used to discover its accelerating expansion (1, 2). They 
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are thought to originate from a thermonuclear explosion of an accreting carbon-oxygen white 
dwarf (WD) near the Chandrasekhar mass limit (Mch ~ 1.4 M ) in a close binary system (3, 
4). Two competing scenarios have been proposed for the progenitor systems: single-degenerate 
(SD) (5, 6) and double-degenerate (DD) models (4, 7). In the former, the mass-donating star 
could be a main-sequence (MS)/subgiant star (8), a red-giant star (RG; 9), or even a helium star 
(10, 11), while it is another WD in the latter scenario (4, 7). Recent results suggest that both 
scenarios are possible (12-18). 

There is increasing evidence for spectral diversity among SNe la. Of particular interest 
are those showing higher expansion velocities as inferred from the blueshifted Si II 615 nm 
feature in optical spectra (19). These fast-expanding SNe la also generally exhibit a steep Si II 
temporal velocity gradient (20). This spectral difference in velocity or velocity evolution of the 
ejecta has been proposed to be a geometric effect of an asymmetric explosion (21, 22). Given 
a common origin for SNe la having different ejecta velocities, they should be found in similar 
stellar environments. This can be tested by examining SN positions in their hosts, the surface 
brightness at these locations, and the properties of their hosts. 

We conducted such an analysis with a well-defined SN sample having 188 SNe la [see sup- 
porting online material (SOM) text SI] from the Lick Observatory Supernova Search (LOSS; 
23). The SN la sample consists of 123 "Branch-normal" (spectra scopically normal) objects 
(24), 30 peculiar ones of the SN 1991bg variety (25), 13 peculiar ones like SN 1991T (26, 27), 
and 7 peculiar ones like SN 2002cx (28), with respective fractions of 65.4%, 16.0%, 6.9%, and 
3.7% (Table SI). There are 15 SNe la (8.0% of all) that cannot be subclassified due to an ab- 
sence of early-time spectra. We concentrate on the Branch-normal SNe la which are thought 
to be relatively uniform. We obtained the main parameters of the host galaxies from two large 
online astronomical databases: the NASA/IPAC Extragalactic Database (NED; 29) and Hyper- 
Leda (30). 
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The location of a SN in its host galaxy can be estimated by the radial distance of the SN 
from the nucleus (Rsn)- Assuming that the galaxies are circular disks and only appear to have 
different major and minor axes due to their inclination, i? SN can be calculated if we know the 
position angle and the axial ratio of each galaxy. The radius of the galaxy (i? ga i) is simply the 
semimajor axis at the 25.0 5-mag arcsec~ 2 isophote. The ratio -RsN/-R g ai is then the fractional 
radial distance of the SN. For SNe la in elliptical galaxies, no tilt correction is applied because 
these galaxies can be regarded as spheroids. The typical host galaxy of our sample (Fig. SI) 
has a major axis of about 1.3'— 1.4' and can be measured with a precision of ~ 0.1'. This results 
in a typical uncertainty of ~ 0.05 in the determination of RsN/R ga ,h 

We measured the velocity of the Si II 615 nm line for 165 SNe la (out of 188) by using the 
published spectral datasets (31-33). We normalized this velocity to the maximum-light value 
with a series of templates of Si II velocity evolution established from well-observed SNe la, 
with a typical uncertainty of 300-400 kms -1 (SOM text SI and Fig. S2). 

It is clear that Branch-normal SNe la with v Si n < 1 2,000 km s _1 (the normal- velocity group, 
NV) span a wide radial distribution, occurring at places from the innermost region to about 2-3 
times the optical radius of the entire galaxy. In contrast, those with v Si n > 12,000 km s _1 (the 
high- velocity group, HV) are rarely found at large galactic radii (Fig. lb). For example, only 3 
out of the 40 HV SNe la are detected in regions with i? SN / R ga \ > 0.7 (2 of which are in elliptical 
galaxies), whereas 14 ± 2 would have been expected at the detection rates of the NV SNe la 
(which are about 34 ± 5% at RsN/R ga \ > 0.7). Such a difference has a statistical significance 
of about 5(7, highlighting the paucity of HV SNe la in outskirts of galaxies. Binning the data 
in velocity space with an interval of about 2000 km s _1 further shows a correlation between the 
ejecta velocity of SNe la and the locations in their host galaxies (see gray quadrangles in Fig. 
la). 

To better understand such a birthplace vs. ejecta-velocity relation for SNe la, it is important 
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to know how the Si II velocity itself is distributed. Fig. lc shows that most of the sample 
of Branch-normal SNe la clusters at velocities between 10,000 and 12,000 km s -1 , with a tail 
extending up to ~ 16,000 km s _1 . Such a velocity distribution can be fit by a double-Gaussian 
model. One component, with a stronger and narrower peak at 10,800km s _1 , is responsible for 
the NV group, whereas the other component, with a weaker and broader peak at 13,000 km s _1 , 
accounts for the HV group. Adopting a velocity cut of t> S i n = 12,000km s _1 to divide these 
two groups puts 40 SNe la in the HV group and 83 in the NV group. Accordingly, we estimate 
the fraction of the HV population to be 1/3 of the Branch-normal sample and 1/5-1/4 of the 
entire SN la sample. Note that the Si II velocity distinction is not sharp between these two 
groups, so blending could occur to some extent. Nevertheless, this blending is small at larger 
velocities, and it will not affect the result that the SNe la with higher velocities tend to occur 
nearer the galaxy centers. 

HV SNe la show the highest central concentration among the samples, with 90% occurring 
in regions within R S N/R ga \ < 0.7; this fraction is 66%, 77%, and 89% for NV SNe la, SNe II, 
and SNe Ibc, respectively (Fig. 2). The contrast is even more significant when examining the 
distribution toward the galaxy center. We caution, however, that a higher fraction of CC SNe 
may be missing in the central regions of galaxies during LOSS because of their lower luminosity 
relative to SNe la. 

A Kolmogoroff-Smirnoff (K-S) test finds a probability of 0.5% that HV and NV SN la 
groups have a similar radial distribution in galaxies. This probability further decreases to 0. 1% if 
one increases the velocity cut dividing these two groups from 12,000 km s _1 to 13,000 km s _1 . 
Possible selection effects in the radial distribution have been explored and none can account 
for such a significant discrepancy between HV and NV SN la groups (SOM text S2 and Figs. 
S3). Thus HV SNe la have higher metallicities and their progenitors are therefore less likely 
to be from the halo population that consists of old, metal-poor stars located far away from the 
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galactic center (34). 

Because the radial distances give only a rough estimate of the properties of SN progenitors, 
more sophisticated methods are necessary to provide additional constraints. A simple statistic 
of the "fractional flux" allows a measurement of how SNe are distributed within their hosts (35). 
This can be achieved by measuring in the host-galaxy images the fraction of total galaxy light 
contained in pixels fainter than or equal to the light in the pixel at the location of the SN. The 
"fractional flux" was obtained with the SDSS u'g'r'-band images (36) for 64 SNe la (39 NV + 
25 HV), 102 SNe II, and 39 SNe Ibc of the entire SN sample (SOM text S3). 

The surface brightness of SN II locations has an approximately linear distribution, roughly 
tracing the distribution of the light in their hosts (Fig. 3). SNe Ibc seem to be more concentrated 
in brighter regions of the host galaxies, consistent with the knowledge that they arise in larger 
star- forming regions that produce more-massive stars. Of the SN la sample, the locations of 
the HV group and the NV group track their hosts' light differently, with a very low probability 
(P = 0.04 in u', 0.07 in g', and 0.08 in r') that they come from the same stellar populations. 
Figure 3 shows that the NV SN la locations are apparently fainter than those of CC SNe; the 
HV SN la distribution, on the other hand, is similar to the CC SN distribution. 

The light radii and luminosities of galaxies in the HV sample and the NV sample also differ 
significantly (Fig. 4). A two-dimensional K-S test gives a probability of 0.5% that they come 
from the same population. The mean-light radius of the hosts estimated for these two groups is 
20.67 ± 0.83 kpc for HV and 16.61 ± 0.67kpc for NV. The mean absolute if-band magnitudes 
are -24.53±0. 13 mag (HV) and -24. 17±0. 10 mag (NV), with the HV hosts being, on average, 
brighter by about 40%. In general, the HV SNe la tend to occur in larger and more luminous 
hosts, and the fraction found in galaxies with R ga \ < 15 kpc is very low, 15% (versus 46% for the 
NV counterparts). This difference is not due to an observational bias because the host galaxies 
of these two groups do not show significantly different distributions in either morphologies or 
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redshifts (see SOM text S4 and Figs. S4 and S5). 

Our analysis thus reveals significant differences in the progenitor environments of SNe la 
having different Si II velocities. The HV SNe la are much more concentrated in regions close 
to the galaxy center and in bright regions of their host galaxies, and they also tend to reside 
in larger and more luminous hosts relative to the NV group. It is generally accepted that all 
galaxies (on average) have metallicities that systematically decrease outward from their galactic 
centers, and that their global metallicities increase with galaxy size (39). Higher metallicity is 
therefore expected for the HV SN la progenitor population. Meanwhile, the fact that the surface 
brightness at HV SN la locations roughly traces the light of their host galaxies, as with CC SNe, 
suggests that the progenitor populations are relatively young. Thus, the HV SN la population 
may have a younger and more metal-rich progenitor system than the NV SN la population, and 
a larger initial MS mass of the exploding WD may be expected for the former due to the shorter 
evolutionary time needed before explosion. 

Calculations of stellar evolution show that both stellar mass and metallicity have significant 
effects on the nature of C-0 WDs that may become progenitors of SNe la. The maximum 
MS mass to form C-0 WDs is found to increase significantly toward higher metallicity (40); 
stars with MS masses up to 8-9 M produce massive WDs (~ 1.1 M ) that can reach 1.4 M Q 
in a shorter time via mass transfer from a companion star. Thus, the HV SN la group might 
represent the young-population SNe la, corresponding to the "prompt" component with short 
delay times (41, 42), while the NV group may belong to the "older" component with long delay 
times. We note that 5 out of the 40 (12.5%) SNe la in the HV sample are in elliptical galaxies 
that are luminous (see Table SI and Fig. S4) and massive. However, among these, some (e.g., 
SN 2002dj and SN 2000B) show dust structure and molecular gas 43, 44), suggestive of recent 
or ongoing star formation in them; they do not contradict the conclusion that HV SNe la likely 
arise from young stellar populations. 
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Having young and metal-rich progenitors may explain the numerous recent detections of 
circumstellar medium (CSM) signatures among HV SNe la. Time-variable absorption features 
of the Na I doublet (Dl 589.6 nm and D2 589.0 nm), likely suggestive of changes in CSM 
ionization due to a variable SN radiation field, have been reported for a few SNe la such as SNe 
2006X, 20071e, and 1999cl (12, 45, 46). A common feature of these SNe la is that they belong 
to the HV subclass (see Table SI). Additionally, the velocity structure of the line-of-sight Na I 
lines provides another possible diagnostic of CSM around SNe. A trend of blueshifts was found 
among SNe la based on a study of Na I absorption lines of a larger SN sample (13). It was 
noticed, however, that the SNe la with a blueshifted absorption feature generally have higher 
Si II velocities (47). This evidence is consistent with their systematically redder color around 
maximum light (perhaps due to additional CSM absorption) relative to the NV population (19), 
and it can be understood in terms of an empirical metallicity dependence of mass outflow (48). 
At higher metallicity, stars lose more mass and produce stronger outflows than their lower- 
metallicity counterparts. Thus, detection of abundant outflows in the vicinity of the HV SNe la 
is not unexpected. 

The higher Si II velocity seen in HV SNe la may be related in part to an increase in stellar 
metallicity. As the metallicity increases in the C+O layer of the exploding WDs, the line- 
forming region moves toward shallower parts of the atmosphere because of an increased line 
opacity, leading to larger line velocities (49). In addition, the high- velocity Si II layers are 
perhaps formed due to the density increase caused by interaction between the SN ejecta and the 
material around the exploding WD, which could be an accretion disk or a filled Roche lobe (50). 
The angular variations in observing such an interacting system may explain the variation of the 
polarization across the Si II line (51) and its correlation with the ejecta velocity for HV SNe la 
(22). This is consistent with the result that the observed differences in some SNe la might be 
due to a projection effect (21). However, these data and analyses are insensitive to a separate 



7 



population of NV SNe la that is intrinsically different from HV SNe la. 

The dependence of SN la ejecta velocity on progenitor environment could be relevant when 
using SNe la as cosmological yardsticks, because the HV and NV populations have different 
colors around maximum light (19) and their ratio may change with redshift. The observed 
relative fraction of the HV and NV population might become smaller at great distances due to 
a decrease in the HV SN la rate in low-metallicity environments and the increased difficulty of 
spectroscopically classifying SNe in the central regions of distant galaxies. 
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Fig. 1: Relations between the Si II velocity of SNe la and the birth location in their host 
galaxies, a, The Si II velocity obtained around £>-band maximum light (v^ n , ordinate) as com- 
pared to the fractional radial distance in the host galaxy (RsN/Rg&i, abscissa) for 165 SNe la. 
The Branch-normal SNe la with n > 12,000 km s" 1 (HV group), those with ugj n < 12,000 
km s" 1 (NV group), SN 1991T-like SNe, and SN 1991bg-like SNe are shown by red triangles, 
blue dots, purple stars, and dark-gray squares, respectively. The SNe la in spiral and ellip- 
tical/lenticular galaxies are represented with filled and open symbols, respectively. The gray 
quadrangles show the radial distances averaged in binned velocity space, which are 0.63 ± 0.08 
in 9,000-12,000 km s \ 0.42 ± 0.06 in 12,000-14,000 km s -1 , and 0.34 ± 0.06 in ~ 14,000- 
16,000 km s _1 , respectively. The horizontal and vertical dashed lines mark the place with vsi n 
= 12,000 km s _1 and with RsN/Rg&i = 0.7. b, The number distribution of the fractional radial 
distance. The red and blue areas are for the HV and NV groups of SNe la. The purple and 
dark-gray areas are for the SN 1991T-like and SN 1991bg-like SNe la. c, The number dis- 
tribution of near-maximum-light Si II velocity. A double Gaussian function is used to fit the 
distribution of 123 Branch-normal (HV + NV) SNe la. Red and blue curves are for the high- 
velocity and normal- velocity components, with respective peaks centered at 13,000 km s^ 1 and 
10,800 km s _1 . The black curve represents the combined result of these two components. 
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Fig. 2: A plot of the cumulative fraction of our SN samples (HV SNe la, NV SNe la, SNe II, and 
SNe Ibc). The gray solid curve represents the distribution of SNe la with t> Si n > 13,000 km s _1 . 
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Fig. 3: Histogram distribution of the fractional flux of the host-galaxy light at the location of 
SN explosions in the vl (left panel), g' (middle panel), and r' (right panel) bands. The diagonal 
black lines represent the case that the SN progenitors follow exactly the distribution of galaxy 
light. The HV SN la population differs significantly from the NV SN la population at the SN 
location (with a respective probability of 3.8% in u', 6.8% in g' , and 8.3% in r' that the SNe 
come from the same radial distribution), but seems to have a distribution similar to those of 
SNe II and SNe Ibc. 
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Fig. 4: A comparison of the physical sizes (£>-band light radius at 25 mag arcsec -2 isophote) 
and absolute f^-band magnitude distributions of the SN hosts. These two parameters have been 
derived by adopting Q m = 0.27, £!a = 0.73, and H = 73kms~ 1 Mpc -1 (37). Foreground 
Galactic absorption corrections (38) have been applied to the absolute magnitudes. In the main 
panel, the high-velocity SN la hosts are represented as red triangles, and the normal-velocity 
SN la hosts as blue dots. The absolute magnitudes of the hosts are shown on the abscissa, 
and the lengths of the semimajor axes of the hosts on the ordinate. The plot is then projected 
onto the two side panels where a histogram is displayed for each host population in each of the 
dimensions, absolute magnitude and semimajor axis. 
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Supporting Online Material 



SI. The Supernova Sample 

The SN sample used in our study is the "season-optimal" subsample from the Lick Ob- 
servatory Supernova Search (LOSS), found over the interval 199872008 (23). This subsample 
contains 499 SNe that were discovered only during (not prior to) the active monitoring pe- 
riod of their host galaxies, with those in small early-type and edge-on spiral galaxies further 
excluded (see Table 4 of Ref. 23). Fig. SI shows the typical host galaxies of the "season- 
optimal" subsample used in our study. Among the 499 SNe, 188 are SNe la (37.7%), 214 are 
SNe II (43.1%), 88 are SNe Ibc (17.6%), and 9 have no known spectroscopic classification 
(1.8%). These fractions of different SN types are fully consistent with those from the full "op- 
timal" sample (726 in total, of which 284 are SNe la) used for the final rate calculations for 
nearby SNe (52), arguing against a possible bias in the sample selection. The main advantage 
of adopting the smaller "season-optimal" subsample is that the discovered SNe are generally at 
a relatively young phase, allowing accurate, timely spectroscopic classifications of them. This 
is particularly important for the identification of spectral diversity among SNe la. 

The ejecta velocity can be measured from the absorption minimum of the blueshifted Si II 
635.5 nm line, which is one of the strongest features in early-time optical/near-infrared spectra 
of SNe la and suffers from little blending with other features. In order to measure an accurate 
expansion velocity from the Si II lines, we smoothed the observed spectrum with a locally 
weighted scatter-plot smoothing method (53) to avoid taking a dip of a noise spike in the data as 
a local absorption minimum. As shown in the inset of Fig. S2, the wavelength of the minimum 
in the smoothed spectrum ( min) is used to calculate the expansion velocity using the special- 
relativistic formula. We further compare these velocities obtained within about one week from 
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Fig. SI: Typical galaxies in the sample. Left: SDSS g'— band image of the host galaxy MCG- 
1-53-6 of SN 2002df, a member of the NV population of SNe la. Right: SDSS g' -band image 
of the host galaxy NGC 3978 of SN 2003cq, a member of the HV population of SNe la. Both 
images have a size of 2.8' x 2.5'. 

maximum light to templates of Si II velocity evolution to determine the corresponding values 
at maximum light, v° Si n (see Fig. S2). These comparison templates are established from some 
well-observed objects with published spectral data (19, 31-33, and references therein). The 
values of v° Si n used in this study are listed in Table SI, with a typical uncertainty of about 
300-400 km s -1 , depending on the number, epoch, and quality of the available spectra. 

Following the classification scheme proposed by Wang et al. (19), we divide the 188 SNe la 
into five groups: high velocity (HV), normal velocity (NV), overluminous SN 1991T-like (91T), 
subluminous SN 1991bg-like (91bg), and the extremely peculiar SN 2002cx-like (02cx) SNe 
la, with respective relative fractions of 21.2%, 44.1%, 6.9%, 16.0%, and 3.7%. (8.1% of SNe la 
could not be classified due to the absence of early-time spectra.) The HVNV ratio, derived with 
the "season-optimal" subsample, may suffer some effects from the excluded sample. However, 
these effects cannot be directly addressed because we do not have early-time spectra of most of 
the excluded SNe and we were thus not able to subclassify them. The spectral data used for the 
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Fig. S2: An example of the Si II velocity measurement and normalization used in this study. 
The absorption feature shown in the inset window is Si II 635.5 nm from a spectrum of SN 
2006cp at t = — 6 day from maximum light. The observed spectrum is the solid black curve, 
the locally weighted scatterplot smoothed spectrum is the solid red curve, and the absorption 
minimum is indicated by the position indicated by the down arrow. The measured velocities 
(represented by the black filled dots) are then normalized to the value at maximum light (indi- 
cated by the red filled dot) by comparing with the template curves of the well-observed SNe la. 
The blue curve shows the mean evolution of NV SNe la, and the gray region represents the la 
uncertainty; the red curves are part of the HV sample; the dashed-plus and dashed-cross lines 
illustrate the evolution of the mean velocity for SN 1991T-like and SN 1991bg-like events, re- 
spectively. These templates are established with the published spectral data from the literature 
(2-4). The black dashed curves represent the interpolated HV templates from the neighboring 
observed ones, with an assumption of proportional change. 
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subclassifications in this study are based primarily on the datasets published by the Center for 
Astrophysics Supernova Program (31, 33) and the Berkeley Supernova Program (32). For some 
objects the spectral parameters are taken from the IAU Circulars. 

Table SI lists the observed parameters as well as the classifications of the 188 SNe la. The 
meaning of each column is as follows: Column (1), the SN name; Column (2), the host-galaxy 
name; Column (3), the recession velocity (v) in the reference frame of the cosmic microwave 
background (CMB) when v > 3000 km s~\ or for a self-consistent Virgo-centric infall vector 
of 220 km s^ 1 when v < 3,000 km s~\ using the procedure given at the NASA Extragalactic 
Database (NED); Column (4), the host-galaxy numerical Hubble type, or "T" type (54); Column 
(5), the host-galaxy absolute K-band magnitude taken from The 2MASS Redshift Survey (55), 
dereddened by the Galactic reddening using the full-sky maps of dust infrared emission (38) 
and the Cardelli et al. reddening law with Ry = 3.1 (56); Column (6), the host-galaxy physical 
radius at the 25.0 fi-mag arcsec~ 2 isophote; Column (7), the deprojected galactocentric distance 
of the SN in its host galaxy, in units of the galaxy radius R ga i; Column (8), the velocity of Si 
II line measured around maximum light, and its uncertainty; Columns (9)-(ll), the u'— , g'— , 
and r' —band fractional fluxes, measured with the SDSS (36) images (57), representing the sum 
of counts registered in all pixels with fewer counts than measured at the SN location divided by 
the sum of all counts associated with the galaxy (35, 58); and Column (12), the SN la subtype. 

S2. Some Selection Effects in the Determination of the Radial Distance 

In calculating the SN relative distance, we assumed that galaxies are circular disks and that 
they exhibit different major and minor axes only because of different inclinations. This assump- 
tion can be used to derive the deprojected radial distance of the SNe in spiral galaxies with an 
intermediate inclination angle (e.g., < 70°. For galaxies with large inclination angles, such a tilt 
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correction may cause large uncertainties in an attempt to derive the deprojected distance. As an 
alternative, we can also make an assumption of spheroids for all of the galaxies. Based on the 
projected radial distance, we reanalyzed the data and found that the HV SN la sample is still 
highly concentrated near the galaxy center relative to the NV sample, with a K-S probability of 
0.3% that they come from the same parent population. 

We further considered possible selection effects that may bias the radial distribution of the 
HV and NV SNe la, which can arise from differences in SN luminosity and size of the host 
galaxy. For example, brighter SNe are more easily detectable in the central regions of galaxies 
than fainter SNe. In Figure S3, we compared the distribution of the V-band magnitudes M^ ax , 
corrected for Galactic reddening, for the HV and NV SN la samples. The mean value and the 
distribution of M^ mx are found to be comparable for these two samples. This indicates that the 
selection effect due to the different luminosity distributions of HV and NV SNe la should not 
have a significant bias on the radial distribution in their host galaxies. 

According to the analysis of the LOSS SN sample (23), there is a higher probability that the 
SNe tend to be missed in the central regions of smaller galaxies relative to that of larger galaxies. 
We further examined this possible selection effect by excluding the smaller galaxies from our 
analysis. Restricting the sample used in the analysis to the SNe la with major axis of the host 
galaxy greater than 1.0' (a total of 96) and 1.3' (a total of 56), we found that the differences in 
radial distributions between the HV and NV subsamples are still statistically significant, with 
K-S probabilities of 0.4% and 5.7% that they come from the same population. 

The above analysis indicates that the concentration of HV SNe la in the central regions of 
galaxies cannot be explained with any known selection effects, so it seems to be a property 
intrinsic to these objects. 

S3. Measurements of the Fractional Flux at the SN Location 
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Fig. S3: Left panel: The V-band absolute magnitudes, corrected for Galactic reddening, plotted 
as a function of the Si II velocity around maximum light. Right panel: Histograms of V-band 
absolute magnitudes for the HV and NV subclasses of SNe la. Two vertical lines in the panel 
mark the respective mean values of the absolute magnitudes for the NV (solid) and HV (dashed) 
SNe la, which are 18.71±0.08 mag and 18.64±0.12 mag. A K-S test rejects the hypothesis with 
a probability of 99.6% that these two SN la populations have different peak luminosities. 
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The fractional flux is measured in host-galaxy images as the fraction of total host light in 
pixels fainter than or equal to the light in the pixel at the location of the SN (35). Were the 
SNe of different types to exactly track the host-galaxy light, their histograms would follow the 
diagonal line. This parameter is independent of galaxy morphology and provides a better way 
to quantify the correlation between SNe and the light of their hosts. The host-galaxy images in 
our analysis were taken from the SDSS DR 6 (56), which includes 54-s exposures in the Sloan 
filter set. To avoid possible contamination by residual light from the SNe, we excluded those 
images in which the SDSS observation of the host was made during the period from 1 month 
before to 12 months after the SN detection. The SDSS u'g'r'— band images were used for the 
measurements; they are available for 25 HV SNe la, 39 NV SNe la, 39 SNe Ibc, and 102 SNe 
II. 

Following Ref. (57), we identify the pixels with host-galaxy light by selecting contigu- 
ous pixels having signal-to-noise ratios greater than 1. Before that, the foreground stars were 
masked out and the enclosed pixel values were replaced with the mean of pixel values around 
the mask perimeter. From the distribution of pixel values of the galaxy light in the SDSS u'g'r' 
bands and the corresponding values at the SN location, we calculate the fractional fluxes. All 
of the results are listed in columns (10)-(12) of Table SI. Note that the u'-band fractional flux 
may suffer larger uncertainties because the SDSS images are not deep enough in this band to 
allow accurate measurements. 

S4. Redshift and Morphology of the Host Galaxies 

Figure S4 shows the redshift distribution of the host galaxies. The mean redshift distribution 
is estimated as z = 0.0240 for the HV SNe la and 0.0217 for the NV SNe la; the HV SN la hosts 
are, on overage, at slightly larger distances. Such a small difference cannot account for the 
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Fig. S4: A plot of the Si II velocities of SNe la near maximum light and the redshift of their 
host galaxies. In the upper panel, the HV SN la hosts are represented as red triangles, the NV 
SN la hosts as blue dots, the 9 IT- like hosts as purple stars, and the 91bg-like hosts as dark- 
yellow squares. The lower panel shows a histogram for the redshift distribution of the hosts of 
different types. The mean redshift for the HV SN la sample is 0.024, comparable to that of the 
NV sample (mean value 0.022).. 
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Fig. S5: A plot of the Si II velocities of SNe la near maximum light and the numerical Hubble 
types of their host galaxies. The symbols are the same as in Fig. S4. The plot is then projected 
onto the lower panel, where a histogram distribution is displayed for each Hubble type of the 
host populations. The K-S test gives a probability of 70% that the HV hosts and NV hosts have 
the same morphological distribution. 
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paucity of the HV SNe la in smaller galaxies, as the mean absolute B-band magnitude and the 
major axis of the host galaxies are found to increase by only ~0.08 mag and ~0.8 kpc from z 
= 0.0217 to 0.0240. The Hubble type ("T" type) of the host galaxies of the HV and NV SN la 
samples is also compared in Figure S5. A K-S test gives a probability of 70% that these two 
samples have a similar morphology distribution, indicating that the differences in light radii and 
luminosity of galaxies seen in Figure 4 are not due to an observational bias. 
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Table SI: Relevant Parameters of the SN la Sample 



SN 


Galaxy 


V 3 k 

(km s" 1 ) 


Galaxy 
Type 


M K 
(mag.) 


Rgal 

(kpc) 


RsN^Rgal V°Si n 

(10 4 kms _1 ) 


fu 




f r 


SN 
Type* 


1998de 


NGC252 


4619 


-1.2 


-24.98(14) 


14.43(125) 


1.54(10) 


1.13(02) 








labg 


1998dk 


UGC 139 


3609 


5.1 


-22.45(19) 


15.56(138) 


0.13(01) 


1.27(03) 


0.88 


0.88 


0.83 


IaH 


1998ef 


UGC646 


5019 


2.9 


-23.78(14) 


15.46(187) 


0.15(02) 


1.33(02) 








IaH 


1998es 


NGC632 


2868 


-1.5 


-22.89(23) 


8.86(109) 


0.23(02) 


1.18(02) 


0.64 


0.48 


0.45 


IaT 


1999aa 


NGC 2595 


4572 


4.9 


-24.34(15) 


28.78(214) 


0.33(01) 


1.04(02) 


0.37 


0.60 


0.52 


IaT 


1999ac 


NGC 6063 


3098 


5.9 


-22.60(22) 


10.49(119) 


0.83(06) 


1.00(02) 





0.07 


0.07 


IaT 


1999bh 


NGC 3435 


5295 


3.0 


-23.58(13) 


18.00(144) 


0.25(02) 


0.88(04) 


0.62 


0.62 


0.54 


labg 


1999by 


NGC 2841 


841 


3.0 


-24.68(23) 


13.02(486) 


0.64(03) 


1.01(02) 


0.50 


0.93 


0.87 


labg 


1999cl 


NGC 4501 


1217 


3.4 


-24.87(55) 


17.17(414) 


0.27(02) 


1.26(02) 


0.44 


0.51 


0.48 


IaH 


1999cp 


NGC 5468 


3112 


6.0 


-22.76(22) 


15.76(161) 


0.75(02) 


1.08(04) 








IaN 


1999da 


NGC 6411 


3748 


-4.8 


-24.45(18) 


16.00(235) 


1.11(13) 


1.07(03) 


0.13 


0.10 


0.02 


labg 


1999dg 


UGC 9758 


6782 


-2.0 


-24.65(10) 


16.21(147) 


0.20(02) 




0.81 


0.41 


0.57 


la 


1999dk 


UGC 1087 


4181 


5.3 


-22.71(18) 


11.67(110) 


0.68(04) 


1.26(06) 





0.30 


0.33 


IaH 


1999do 


MCG+05-54-3 


6238 


4.7 


-23.36(12) 


5.42(106) 


0.36(07) 










la 


1999dq 


NGC 976 


4060 


4.1 


-24.65(16) 


11.97(097) 


0.19(01) 


1.07(02) 


0.85 


0.72 


0.64 


IaT 


1999ej 


NGC 495 


3831 


0.1 


-23.66(17) 


9.40(093) 


0.82(05) 


1.09(10) 


0.28 


0.22 


0.22 


IaN 


1999ek 


UGC 3329 


5278 


4.1 


-24.71(13) 


17.66(166) 


0.63(21) 


1.07(02) 








IaN 


2000A 


MCG+01-59-8 


12891 


7.2 


-24.80(07) 


27.74(370) 


0.16(02) 


0.97(05) 








IaN 


2000B 


NGC 2320 


6020 


-4.8 


-25.76(11) 


16.38(291) 


0.56(09) 


1.45(07) 








IaH 


2008hv 


NGC 2765 


4074 


-1.9 


-24.53(16) 


15.47(136) 


0.89(09) 


1.10(10) 





0.26 


0.24 


IaN 



Note: The uncertainties shown in the parentheses are in units of 0.01 mag for M K , 0.01 
kpc for Rg a i, 0.01 for relative distance, and 100 km s -1 for Si II velocity. 

*'Iabg' = 91bg-like SNe la; 'IaT' = 91T-like SNe la; 'IaN' = NV SNe la; 'IaH' = HV SNe 
la, 'IaPec' = 02cx-like SNe la; 'la' = unclassified SNe la. 

The full data of Table SI is available as a separate text file. 
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